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By applying methods ol integrable quantum field theory 
to the Kondo problem, we develop a systematic perturbation 
expansion near the IR (strong coupling) fixed point. This 
requires the knowledge of an infinity of irrelevant operators 
and their couplings, which we all determine exactly. A low 
temperature expansion (ie all the corrections to Fermi liquid 
theory) of the resistivity then follows, extending for instance 
the well known Nozieres T 2 result in the exactly screened case 
to arbitrary order. The example of the ordinary Kondo model 
is worked out in details: we determine p up to order T 6 , and 
compare the result with available numerical data. 

Theoretical progress in the understanding of the 
Kondo model over the last twenty years has been remark- 
able (see jlj for a review). Most of the interest has fo- 
cussed on the strong coupling fixed point, where non per- 
turbative techniques like the numerical renormalization 
group, the Bethe ansatz || and conformal field theory 
H have provided many powerful results. Yet, the ba- 
sic quantity that gave rise to the discovery of the Kondo 
physics itself, the resistivity p, has remained strangely 
resilient to analytical understanding : although its be- 
haviour is well understood, based on various numerical 
renormalization or approximation methods, essentially 
no exact results have been obtained for it - for instance, 
it cannot be computed by the Bethe ansatz, in contrast 
with the specific heat, or the magnetic susceptibility. A 
notable exception to this state of affairs is the result of 
Nozieres Q for the ordinary Kondo problem, which gives 
the exact first order correction to the resistivity near the 
strong coupling (IR) fixed point (see below). This result 
has since been generalized to the case of several channels 
and higher spins |Q. 

It is important here to emphasize that studying the 
vicinity of the strong coupling fixed point is a subtle 
challenge. In the case of the weak coupling (UV) fixed 
point, one at least knows the effective hamiltonian, up 
to irrelevant terms, so the universal properties at high 
temperature can in principle be accessed perturbatively. 
As is well known, this perturbation breaks down at some 
point. A natural idea to study the low temperature prop- 
erties is thus to perturb around the other, strong coupling 
fixed point. In that case however, all the operators are 
irrelevant. The leading irrelevant operator is usually the 
only one kept in the computations: it is this operator 
that gives rise to the Fermi liquid or Nozieres result in 
the exactly screened case, or the results of Affleck and 
Ludwig in the overscreened case ||. To go beyond first 
order, one however needs to know the next irrelevant op- 



erators, together with their coupling constants. This is 
a difficult problem, that cannot be answered by qualita- 
tive arguments. As we will see below, even the form of 
the operators is not obvious, and requires a complicated 
adjustment. The only way to determine precisely the 
hamiltonian near the IR fixed point is to use, somehow, 
the fact that it is the "large coupling limit" of the known 
hamiltonian near the UV fixed point. How to do that in 
practice however was a largely open problem up to now. 
A few years ago, Hewson |6| introduced a "renormalized 
perturbation theory" that should in principle answer the 
question : in practice however, the method seems very 
hard to implement, and we are not aware of any quanti- 
tative correction to the Fermi liquid behaviour obtained 
in that way. 

The main result of this paper is the exact determina- 
tion of the hamiltonian near the IR fixed point to all 
orders. It is then a purely technical matter to determine 
the corrections to the resistivity. We carry them out to 
order T e in the ordinary Kondo problem, and compare 
with existing numerical data. 

Let us now define the problem in more details. We 
follow closely the notations of || . We take as the starting 
hamiltonian 
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with fi the spin label. Assuming that the impurities are 
dilute enough so we can restrict to the case of only one 
impurity, we decompose as usual the fermionic field on 
spherical harmonics. Only the s-wave interacts with the 
impurity and, at low energies, its field operator takes the 
form 
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with the proper linearisation of the dispersion relation 
in that limit, ie e p ~ vp(p — pp) = vpp' . The resistiv- 
ity can be obtained by using the Kubo formula for the 
conductivity, following ]Bj 
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with the single particle lifetime defined by 1/t = 
— 2Im£ ii . the retarded self-energy, is obtained from 
the three-dimensional Green's function which in turns 
follows from the one-dimensional s-wave contributions || 
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Here the superscript means that the Green's function 
is evaluated with respect to the non-interacting UV fixed 
point where i/>_r(0) = iPl{0) as a consequence of the 
spherical harmonic decomposition. Eq (^) can be rewrit- 
ten as 



G(ujM,ri,f 2 ) - G°(uj M ,r 1 - r 2 ) = 
G° (uj m , n )T{cjm )G° {ojm , -r 2 ). 
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One can then sum over the dilute array of random im- 
purities, restoring translation invariance. The summa- 
tion over multi-impurity terms leads to the self energy 
E(wjk) = riiT(ujM), to first order in rii, the impurity 
density. As usual, impurity-impurity interactions are ne- 
glected, and the following results should be reliable for 
dilute impurities (see ||, for more discussion on that 
point). The retarded self-energy is found by analytical 
continuation, iujm —> Lu+iO. It leads to the single particle 
lifetime, and thus to the resistivity using (||). 

It is thus clear, as stressed in ||, that the RL and LR 
components of the ID Green's functions are the only nec- 
essary ingredients to compute the resistivity, and we will 
concentrate on these in the following. At the IR fixed 
point, the spin is completely screened and the interac- 
tion with the impurity amounts to a phase shift, which 
in the one-channel case is simply encoded in the bound- 
ary conditions V'-r(O) = — V'i(O): this leads to the result 
T, R (uj) = — i^ij where v is the density of states per spin. 

To proceed, we bosonize the one dimensional theory, 
setting tpB./L.fj, — exp (iiV^t^) (our normalization for 
the fermions follows || anddoes not contain the usual 
factors of 2ir). The spin and charge degrees of freedom 
then completely decouple from each other and, using spin 
and charge fields 

<t>s = —?=(^T - 0J.)> <t>c = + 0i)> ( 6 ) 

the interaction only involves the spin field, with hamilto- 
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where right and left movers have been combined in one 
field, and U* = \ /° M dx[{d x <j> s ) 2 + Il 2 s ] is a free bo- 
son hamiltonian. To write (0) we have performed an 
additional canonical transformation to get rid of the S z 
terms. In the usual approximation where the band struc- 
ture does not mix up spin and charge, the RG flow is 
completely controlled by the spin sector. 

The hamiltonian (^) defines an integrable boundary 
quantum field theory [pi. Although much progress has 



been made in computing physical properties by using in- 
tegrability |^], including some correlators ||, the resis- 
tivity so far has not been accessible. The two reasons for 
that are: the divergence of the form factors expansion 
when the perturbing operator has dimension 1 as in (|?|), 
and the fact that the resistivity is really a three dimen- 
sional quantity, following in a very non trivial way only 
from the one dimensional computations. 

Our strategy here will be to use integrability to deter- 
mine the exact form of the IR hamiltonian. Once this is 
obtained, we will be able to do perturbation theory near 
the IR fixed point: short of a Bethe ansatz determination 
of the resistivity, this seems the most powerful approach 
at the present time. 

We write 

oo 

H m = H* + 5H = H* + Y, u 2 k + iO 2k+2 {0), (8) 

fc=0 

where the coupling constants u and operators O have to 
be determined. In (||), Tt* is still a fre boson hamiltonian. 
However, due to the modified boundary condition ipL = 
—ipR at the IR fixed point, we have <^l(0) = 4>R{0)+i^/j- 
It is fair to say at this stage that the determination 
of (|J) follows implicitely from the works of Bazhanov 
et al. fllOj . We present here a much more direct, and 
generalizable, argument. The main idea is the following: 
we first restrict to T — 0, and we switch the description 
from "open string channel" as in (Q) (where imaginary 
time is in the y direction) , to the "closed string channel" , 
where imaginary time is now along x. In that case, the 
theory is described by a free boson hamiltonian, and all 
the interactions are encoded in a boundary state \B). 
The latter takes a simple form M because of integrability: 
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Here, |0) is the ground state of the theory in the open 
string channel, Tb is proportional to the Kondo tem- 
perature (we discuss their exact relation below), the Z* 
are creation operators for the kinks and antikinks of the 
massless sine-Gordon theory jll| - that is, the basis diag- 
onalising the interaction. (3 is the rapidity of the mass- 
less kinks, satisfying e = ±p = e@ for right (resp. left) 
movers. We do not need the precise way in which the 
kinks and antikinks behave to explain the argument: it 
suffices to say that their scattering in the bulk is factor- 
ized, and that, in the open channel picture, they bounce 
on the boundary in a factorized way too, with a re- 
flection matrix that is simply related with K. Here, 
K(j3) — — itanh^. To proceed, let us now expand \B) 
around the IR fixed point, ie for Tb large. We get 
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where \B)* is the IR boundary state and the T2k+i are 
the standard conserved quantities of the theory: they act 
diagonally on the multiparticles states with eigenvalues 
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with d the chirality. By comparing with the general 
formula \B) = Pexp[— J dySH]\B}* , where Pexp is the 
y-path ordered exponential, we can identify the hamilto- 
nian. The last step to do so is to reexpress the I2/C+1 in 
terms of (integrals of) local operators of the theory. This 
is easy to do, using the condition that these operators 
must be mutually commuting (how to find the proper 
normalization will be discussed in a subsequent publica- 
tion |rl| ). After a Wick rotation to go back to the open 
string channel, the final result is 
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the irrelevant operators are dictated by the microscopic 
theory, and usually all of the same order, so only the lead- 
ing relevant one remains in the universal scaling regime. 
Here, the vicinity of the IR fixed point is fully determined 
by the constraint that it is the large coupling limit of the 
UV one, and this is what produces the fine tuning. 

It is important to emphasize that our IR hamilto- 
nian makes sense only within a well defined regulariza- 
tion scheme. Since all this is based on the formalism of 
integrable quantum field theories where no cut-off ever 
appears, we have to regularize our integrals by contour 
splitting, as explained below. 

Let us now discuss how to compute the resistivity. The 
first operator perturbing the fixed point is : (d y 4> s ) 2 :, the 
energy momentum tensor. To order ^ for instance, the 
correction to the RL ID Green's function for the up-spin 
electrons reads 

SG ( nl = 7f- / dye^'y / dy' (15) 

J B J -(3/2 J -0/2 
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Here, the basic object is the stress energy tensor T ww — 
— 2ir : (d w 4> s ) : (w — —y + ix); it is the leading irrelevant 
operator near the IR fixed point, a result well known from 
previous work J5|. Up to an overall normalization, there 
is no ambiguity it: it is the only operator of dimension 
2, up to total derivatives. This is hardly true however 
for the next to leading irrelevant operators: they involve 
very particular combinations of powers of derivatives of 
cj> s , and could not have been guessed on general grounds, 
without using the constraint of integrability. The case 
of O2 is an exception, where, after refermionization, one 
finds 

2 oc 2 : VjV'rV'IV'l : ~i ■ ^\d y ^ : -i : ip\d y tpi : (14) 

It is proportionnal to J(0) • J(0) (up to total derivatives), 
with J the spin current fl4|| . 

Notice that the more irrelevant operators in come 
with prefactors that are inverse powers of the Kondo tem- 
perature: due to this rescaling, they cannot be neglected, 
and a perturbation expansion to n th non trivial order 
near the IR fixed point requires the knowledge of the 
terms in up to k = n. This is in contrast with the sit- 
uation near the UV fixed point, where the coefficients of 



where the propagators have to be evaluated with the IR 
fixed point hamiltonian. The integrals in ([l5]) can be 
quickly evaluated by closing the contour in the complex 
plane using the periodicity of the integrand. The result 
is simply 
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Since this term is purely imaginary, it contributes a real 
term to the self energy, and does not affect the resistivity. 

At order O2 is still the only operator to contribute. 

B 

The integrals now present divergences where two of the 
O2 operators come close. These divergences are simply 
regulated by shifting one of the contours slighlty off in 
the imaginary direction. One finds finally 
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in agreement with g. 

We now extend the computation beyond this order. At 
order 1/Tg, it becomes necessary to include, in addition 
to the energy momentum tensor, the higher conserved 
operator O4. In manipulating O4, one has to be very 
careful that it is not a primary field: its correlation func- 
tions on the cylinder (ie at finite temperature) contain 
anomalous terms, that can be computed from general 
conformal field theory arguments (for instance, one finds 
: Tl w := 47T 2 : {d w cj> s ) 4 : -2tt : d w <fi s dl<j> s : -2tt 3 T 2 : 
(d(j) s ) 2 : + 3w 2 q )• ^ s before, this order does not con- 
tribute to the resistivity. At order 2 and O4 are 
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still the only operators that contribute, and we find 
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In these computations, very strong short distance diver- 
gences are of course generally encountered (there are no 
large distance divergences due to the temperature T that 
acts as a cut-off). To regularize them, as for the second 
order case, we slightly move the contours off one another, 
and use the residue theorem to evaluate the resulting in- 
tegrals. The commutativity of all the l2k+i means that 
the short distance expansions of the operators 02fc+2 
have simple poles whose residues are total derivatives: 
this ensures that the final result does not depend on the 
order in which this regularization is implemented. 

We have carried out the computation further to fifth 
and sixth order. The expression for the self energy is too 
bulky to be given here, but, by using the Kubo formula 
as in [pj, it leads to our main result for the resistivity 
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where v is the density of states per spin. The order T 2 al- 
lows us to quickly match normalizations with other works 
- see for instance appendix K in the book by Hewson |l]] 
and we find Tb = Tk the usual Kondo tempera- 

ture, w = 0.41071 the Wilson number (compared to the 
review by Andrei et al. g, T B = 2T ). The orders T 4 
and T 6 in © are new. 

In general theories, one expects IR expansions to be 
only asymptotic. However, for integrable quantum im- 
purity problems, these expansions often turn out to be 
convergent. Assessing the status of the T 2 expansion of 
the resistivity is a bit difficult: based on our numerical 
values, it seems to converge, but with a disappointingly 
small radius of convergence. That the radius is so small 
presumably explains why estimates of Tk from numer- 
ical or experimental data on the resistivity often differ 
significantly from estimates based eg on the susceptibil- 
ity (where convergence is good, as is easily checked on 
the Bethe ansatz results). The comparison of our for- 
mula (^|) with the result of numerical renormalisation 
group methods carried out by Costi et al. Eq] on the An- 
derson model is shown on the following figure. Clearly 
more orders would be needed to go beyond T/T B w .2. 
Pade approximants however turn out to be quite stable, 
and agree well with the data up to T/Tb ~ 1, which is 
well in the crossover regime. We thus believe that our re- 
sults present a rather complete analytical understanding 
of the strong coupling resistivity. 



1 .2 

1 .1 
1 

0.9 
0.8 
S 0.7 
§0.6 
c£0.5 
0.4 
0.3 
0.2 
0.1 



Numerical RG 
IR perturbation 
Pade 





10 




FIG. 1. Comparison of IR perturbation theory with numer- 
ical RG results. 
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